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Using a spin-rotation invariant version of the slave-boson approach we investigate the relative stability and 
band structure of various incommensurate phases in the cuprates. Our findings obtained in the Hubbard 
model with next-nearest neighbor hopping —t'/t ~ 0.15, as appropriate for the La2-2,Sra;Cu04 family, 
support the formation of diagonal (vertical) stripe phases in the doping regime x — 1/16 (a; — 1/8), 
respectively. In contrast, based on the fact that a larger value —t'/t = 0.3 expected for YBa2Cu306+i 
triggers a crossover to the diagonal (1,1) spiral phase at increasing doping, we argue that it might explain 
why the static charge order has been detected in YBa2Cu306+5 only in the highly underdoped regime. 

Copyright line will be provided by the publisher 

The abundance of experimental results and theoretical investigations in recent years have shown that the 
ground state of the high-T^ cuprates might be spatially inhomogeneous [1 1. In particular, both charge and 
spin modulations have been detected in neutron scattering studies performed on Lai g_2,Ndo.4Sra;Cu04. 
Remarkably, incommensurate (IC) Bragg peaks indicate that in the lightly doped regime x < 1/16 the 
antiferromagnetic (AF) domains are separated by diagonal charge domain walls (DWs) |2| but at a higher 
doping level a; = 1/8 the doped holes self-organize into stripes that run vertically across the Cu02 planes 
[31 . On the one hand, the former are on average filled by 1 / \/2 hole per two DWs or even by one hole 
per one atom in a DW as established in La2_a;Sr2;Cu04 (LSCO) around x — 0.02 |4|, which corresponds 
to the so-called ^//ec/ stripes. On the other hand, the latter are characterized by the filling of one hole 
per two atoms along the DWs meaning that they are only half-filled. IC charge order consistent with 
vertical/horizontal stripes has also been reported in highly underdoped YBa2Cu306+5 (YBCO) with S = 
0.35 [5], but despite several attempts no static charge order could have been detected at higher doping. 
Apart from explaining the neutron scattering data, stripe phases also offer a framework for interpreting 
a broad range of other experimental results such as angle-resolved photoemission spectroscopy (ARPES) 
||6| (see below) as well as magnetic excitation spectra in both LSCO and YBCO families that might be 
consistently explained in terms of fluctuating stripes fV.^. 

While the stripe debate continues, an alternative scenario which might account for the IC spin structure 
in LSCO is a deformation of the AF order so as to optimize the hole motion that stabiUzes a spiral phase 
milPl. Moreover, Lindgard flT] has shown that spiral states can also resolve the universality of magnetic 
excitations in the cuprates and provide a competing paradigm with the stripe phase concept. Therefore, 
the main purpose of this paper is to address in a systematic way the competition between two possible 
site-centered (i.e., with a fully suppressed magnetization along the domain walls) stripe ground states: 
half-filled vertical site-centered (HVSC) and filled diagonal site-centered (FDSC), their bond-centered 
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Fig. 1 Effect of the next-nearest neighbor hopping t' on the density of states at doping a; = 1/8 as obtained in the 
extended Hubbard model for the uniform: (a) AF phase; (b) (0,1) spiral phase, and (c) (1,1) spiral phase. Parameters: 
U = 12t, t' = (dashed line), and £' = -Q.3t (solid line). 



counterparts (where DWs are given by ladders with a weak ferromagnetic order on the rungs): half-filled 
vertical bond-centered (HVBC) and filled diagonal bond-centered (FDBC) DWs, as well as between the 
vertical (0,1) and diagonal (1,1) spiral phases. In order to treat noncanted stripe phases and the spiral order 
on equal footing and to implement local electron correlations we employ a rotationally invariant version of 
the slave-boson (SB) approach in spin space fT2l by introducing auxiliary boson operators {e,;, di,pia,Pj} 
which control the actual electronic configuration at each site i. The Hamiltonian may be then written as, 

ij a\r i 

where {z j , } are 2x2 matrices in spin space which depend on the actual configuration of the boson fields, 
the hopping <y is t on the bonds connecting nearest-neighbor sites and t' for next-nearest neighbor sites, 
while U is the on-site Coulomb interaction. Calculations were carried out on square 128x128 clusters 
which became possible by developing an efficient scheme in reciprocal space which makes use of the 
stripe symmetry ifTSll . In our studies, we have chosen U = 12t, which gives the ratio of J/t = 1/3 (with 
the superexchange J ~ At'^/U), being a value representative for LSCO [14]. 

Regarding the excitation spectrum of the Hubbard model, dynamical mean-field theory (DMFT) re- 
vealed that it consists of a lower Hubbard band (LHB), an upper Hubbard band (UHB), and a quasiparticle 
peak that vanishes at the Mott transition 1 1 5 1 . Since the former two are dynamical in nature, they are hardly 
present in mean-field calculations, while the latter is at the heart of the Brinkman-Rice description of the 
Mott transition and of the equivalent SB mean-field approach in the paramagnetic phase. In contrast, in 
stripe phases bands, separated by energies of order U and corresponding to the LHB and UHB, are part of 
the spectrum, on top of midgap bands in which the chemical potential may be located, as discussed below. 

Let us now review the mechanisms responsible for the stability of these various IC phases in the pure 
Hubbard model with t' = at moderate doping x < 0.15. First of all. Fig. [T] shows that a prominent 



Table 1 Local magnetization rrii and double occupancy di as found for the uniform AF state as well as for (0,1) and 
(1,1) spiral phases in the extended Hubbard model with U = 12t for x — 1/16, and x — 1/8. 
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1/16 


-0.3 
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X = 
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1/8 
t'/t = 


-0.3 


phase 


rrii 




nil 


di 


■nil 


di 


nil 


di 


AF 
(0,1) 
(1,1) 


0.803 
0.813 
0.820 


0.0265 
0.0257 
0.0251 


0.799 
0.811 
0.829 


0.0257 
0.0247 
0.0229 


0.586 
0.646 
0.682 


0.0227 
0.0216 
0.0202 


0.595 
0.656 
0.730 


0.0207 
0.0190 
0.0158 
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Table 2 Free energy F/t per site for various phases at x = 1/8 found in the SB approximation in the extended 
Hubbard model with U — 12t for representative values of the next-nearest neighbor hopping: t' — 0, t' — — 0.15f, 
and t' — ~0.3t. The most advantageous energy for each t' is given in bold characters. 



t'/t 


AF 


(0,1) 


(1,1) 


FDSC 


FDBC 


HVSC 


HVBC 


0.00 
-0.15 
-0.30 


-0.5393 
-0.5323 
-0.5341 


-0.5613 
-0.5564 
-0.5594 


-0.5700 
-0.5699 
0.5842 


0.5821 
0.5716 

-0.5655 


-0.5819 
-0.5713 
-0.5651 


-0.5689 
-0.5700 
-0.5749 


-0.5680 
-0.5691 
-0.5740 



feature of the spiral phases is the van Hove singularity below the Fermi energy that leads to a substantial 
energy gain with respect to the AF phase. Remarkably, this energy gain is amplified by doping making the 
spiral phases increasingly competitive. Next, as listed in Table[T] they are characterized by larger staggered 
magnetization resulting into smaller double occupancy and consequently into a better Coulomb energy as 
compared to the AF phase. Finally, one finds that this mechanism is more effective when the spiral twist 
is along the diagonal of the Brillouin zone. In contrast, stripe phases result from a compromise between 
the optimized kinetic energy of doped holes and the superexchange interaction |[T3l . Indeed, propagation 
of doped holes is easy when the magnetization is locally suppressed along the charge DWs, while the 
superexchange is best optimized when the AF domains are close to half-filling. 

Turning now to the extended Hubbard model, the main effect of t' is to shift the van Hove singularity 
towards negative energy 1161 . As a result, perfect nesting of the Fermi surface is broken and a metal- 
to-insulator transition occurs for a finite Uc, even at half-filling 1 17|. Accordingly, antiferromagnetism is 
suppressed at a; = becoming, however, enhanced when moving towards the van Hove point. Similarly, 
the filled stripe phases should not benefit from t' as it strongly frustrates the antiferromagnetism in the 
weakly doped AF domains separating DWs. In contrast, one can conjecture that the spiral phases better 
make use of the van Hove singularity by its systematic shift below the Fermi energy. This, together with 
both larger staggered magnetization and smaller double occupancy, should yield a more significant energy 
gain and make the spiral phases more competitive with increasing —t'/t. Therefore, one may expect a 
critical t' for which the (1,1) spiral phase becomes the ground state. As shown in Table |2] this indeed 
happens for t' ~ ~0.2t at doping x ~ 1/8. In this case, one observes a spectacular transfer of density of 
states from the Fermi energy towards the bottom of the LHB as depicted in Fig.[TJc). In addition, staggered 
magnetization in the (1,1) spiral phase is strongly enhanced as compared to the AF phase supported by 
concomitant reduction of double occupancy (see Table [U. However, Table[3]shows that for small doping 
a; = 1/16, the above mechanisms are not sufficient to stabilize the spiral order even for t' /t — —0.3. 

Remarkably, since the holes are more delocalized over the AF domains which are then less frustrated, 
finite t' also promotes the half-filled vertical stripe phase. In addition, it induces conspicuous changes in 
the band structure that should have measurable consequences in ARPES experiments. First of all, as shown 
in Fig. 121 a) segregation of doped holes into metallic site-centered DWs gives rise to the appearance of two 
additional partially filled bands lying within the Mott-Hubbard gap. Secondly, as a consequence of the 
vertical stripe order, the T — X and T — Y directions with X = (tt, 0) and Y = (0, tt) are inequivalent. In- 
deed, as the fully quenched spin polarization facilitates the propagation of the holes along the site-centered 
DWs, the maximum dispersion is expected along the T — Y direction. In contrast, the midgap bands al- 
most do not disperse along the T — X direction, which indicates that the holes are rather static along the 



Table 3 The same as in Table[2]but for a; = 1/16. FDBC stripe phase is unstable at t' = —0.3t. 



t'/t 


AF 


(0,1) 


(1,1) 


FDSC 


FDBC 


HVSC 


HVBC 


0.00 

-0.15 
-0.30 


-0.4183 
-0.4173 
-0.4194 


-0.4286 
-0.4287 
-0.4319 


-0.4363 
-0.4376 
-0.4444 


0.4561 
0.4508 

-0.4478 


-0.4560 
-0.4507 


-0.4496 
-0.4503 
0.4529 


-0.4492 
-0.4499 
-0.4525 
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Fig. 2 Low energy part of the electronic structure found at a; = 1/8 
for the HVSC stripe phase in the extended Hubbard model with 
U = 12t as well as with: (a) t' = and (b) t' = -0.15t. Thick 
solid line depicts midgap states while dashed line indicates the Fermi 
level. The highest energy band is centered around lj — /i ~ 7f. 




Fig. 3 The same as in Fig. |2{b) 
but for the HVBC stripe phase. The 
two highest energy bands are cen- 
tered around cu — ^ ~ 6t. 



horizontal direction. Therefore, both midgap bands could be approximated by one-dimensional noninter- 
acting dispersion e^* = ±2icosfcy + est. It crosses the Fermi energy at k = (7r,7r/4), (7r/4, 7r/4), 
and at the equivalent points as expected for the stripes separated by four lattice spacings. However, due 
to strong local coiTelations the bandwidth of these states is strongly renormalized and reduced down to 
VFdw /t ~ 1.45. On the one hand, some features depicted in Fig.|2a) such as the flat band near the X point 
crossing the Fermi level at the k — (tt, 7r/4) point and excitations at the S — {tt/2, 7r/2) point at a higher 
binding energy that at the X point, are in a close agreement with the ARPES spectral density for LSCO at 
X — 1/8 [6|. On the other hand, no spectral weight was found at the Fermi energy along the nodal F — M 
direction with M — (tt, tt), especially at the (7r/4, 7r/4) point. Therefore, one has to use more accurate 
methods, such as the DMFT, which include dynamical coiTelations and leads to a better agreement with 
the experimental data, reproducing a pseudogap at a; = 1/8 [T?]. Furthermore, Fig. |2fb) indicates that 
a negative t' ~ — 0.15t, as expected for LSCO, leads to a distinct broadening of the midgap bandwidth 
up to Wuw/t — 1-6. Consequently, it shifts these states to a lower energy which explains the increasing 
stability of the half-filled DWs in spite of the concomitant narrowing of the LHB bandwidth associated 
with the insulating AF background. Finally, similar conclusions might be drawn concerning the effect of 
t' on the band structure of the HVBC stripe phase shown in Fig. [3] In this case, one finds a gap between 
both midgap states due to a finite AF spin polarization along the bond-centered DWs. However, shape of 
the band structure depends on the charge periodicity. Indeed, even period stripes yield a spectral gap at the 
5* point while this gap is found to be absent for odd period ones ifTSl . 

An entirely different mechanism is responsible for the stability of the filled site-centered DWs. It might 
be best understood from the band structure shown for x — 1/16 in Fig. Ufa). Here, each diagonal DW 
induces the formation of two (dispersionless along the X — Y direction) entirely unoccupied midgap 
bands. Thus the special stability of these phases follows from a real gap that opens in the symmetry broken 
state between the highest occupied state of the LHB and the bottom of the midgap band fT9"20'|. Besides, 
the main effect of increasing — i'/t is to narrow the LHB and to push its lowest energy states towards higher 
energy while the midgap states remain intact as depicted in Fig.|4|b). This feature is clearly visible along 
the parallel to the stripes T — M direction and strongly decreases the stability of the FDSC stripe phase 
with respect to the half-filled one (see Tables |2] and O. Unfortunately, the above gap is inconsistent with 
the ARPES data on the lightly doped LSCO. Indeed, experimentally a quasiparticle peak crosses the Fermi 
level in the nodal direction F — M and forms a hole pocket at the S point fT\\. A possible interpretation is 
provided by Fig. |5] showing that finite spectral weight around the S point may arise from the bond-centered 
DWs. Here, finite spin polarization of DWs leads to the hybridization of the midgap bands with the LHB. 
Such feature was also found in the two-band Hubbard model describing stripe phases in the nickelates 1221 . 
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Fig. 4 Low energy part of the electronic structure found at a; — 1/16 Fig. 5 The same as in Fig. |4{b) 

for the FDSC stripe phase in the extended Hubbard model with but for the FDBC stripe phase. The 

U = 12t as well as with: (a) t' = and (b) t' = — 0.15t. The highest highest energy bands (not shown) 

energy bands (not shown) are centered around a; — p ~ 9t. are centered around — ^ ~ 9t. 



In conclusion, we have found that the next-nearest neighbor hopping t' plays an important role in affect- 
ing the relative stability and the band structure of the IC phases. Next, partial filling of the bands within the 
Mott-Hubbard gap, pseudogap that forms at the Fermi energy around the S point, and the flat band near 
the X point crossing the Fermi level at the k = (tt, 7r/4) point are found to be general consequences of 
the half-filled vertical stripe phase in qualitative agreement with the ARPES results established at doping 
X = 1/8. In contrast, finite spectral weight around the S point experimentally observed at a; — 1/16 
doping might be attributed to the FDBC stripe phase also explaining the metallic behavior of the lightly 
doped LSCO. Finally, we would like to emphasize that the strongly enhanced stability of the (1,1) spiral 
phase by finite t' above a; = 1/8 with respect to the above stripe phases is robust and is not affected even 
by optimizing the stripe filling that further lowers the energy of the stripe phases 123 ]. This may explain 
why the static charge order has been detected in YBa2Cu306+5 only in the highly underdoped regime. 
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